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I.  Introduction 


Compressive  sensing  (CS)  has  received  considerable  research  interest  recently.  CS  has  been 
successfully  applied  in  many  applications,  such  as  high-resolution  image  processing,  wireless 
communication,  electronic  signal  capturing  and  biomedical  imaging  systems  [1],  In  the  context 
of  radar,  the  number  of  potential  targets  is  often  limited.  Hence,  exploiting  the  sparsity  of  radar 
signals  in  various  spaces,  CS  has  been  recently  applied  to  multiple-input  multiple-output  (MIMO) 
radar  systems. 

In  principle,  CS  focuses  on  the  role  of  sparsity  in  reducing  the  number  of  measurements  needed 
to  represent  a  finite  dimensional  vector.  Since  recovery  from  compressed  measurements  is  an  NP- 
hard  problem,  many  suboptimal  methods  have  been  proposed  to  approximate  the  sparse  solution. 
These  algorithms  recover  the  true  value  of  the  sparse  vector  when  the  columns  of  sparsity  basis 
matrix  are  incoherent,  called  mutually  incoherency  (MC)  criteria  [7],  [14],  [8],  [9]. 

A  related  criteria,  referred  to  as  restricted  isometric  property  (RIP),  is  also  used  to  guarantee 
the  accuracy  of  sparse  signal  recovery  [2],  [3],  [15],  [11],  [9],  [16].  However,  in  practice,  it  is 
difficult  to  verify  incoherency  or  RIP  criteria.  It  is  well  known  that  a  random  measurement  matrix 
typically  satisfies  the  incoherency  condition  in  most  cases.  Therefore,  in  related  published  works, 
the  elements  of  the  matrix  $  are  randomly  selected  from  a  Gaussian  distribution  or  a  random 
pseudo  random  sequence  }<g)l|  . 

However,  it  is  difficult  to  determine  the  probability  of  perfect  recovery  in  the  context  of  CS  theo¬ 
ry.  Recently,  some  efforts  have  been  made  to  find  conditions  and/or  tight  bounds  on  the  probability 
of  perfect  recovery  (see  [4],  [5],  [6]  and  the  references  therein).  Although  the  nature  of  mutual 
coherence  is  random,  the  existing  CS  methods  consider  MC  or  RIP  as  constants  when  deriving  the 
probability  of  perfect  recovery  [8],  [9],  [15],  [16]. 

In  this  work,  we  provide  a  probabilistic  approach  for  the  problem  of  perfect  sparse  signal  recov¬ 
ery  under  different  incoherency  conditions.  To  the  best  of  our  knowledge,  this  is  the  first  paper  that 
models  both  MC  and  RIP  conditions  as  the  combination  of  numerous  independent  and  dependent 
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random  variables  (RVs).  In  this  contribution,  our  ultimate  goal  is  to  find  the  perfect  reconstruction 
probability  for  CS -based  MIMO  radars.  Our  results  along  with  the  developed  conditions  in  [16] 
and  [9]  can  be  used  to  estimate  the  maximum  number  of  detectable  targets. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  II,  we  first  present  the  signal  model. 
Then,  we  discuss  about  the  required  conditions  for  the  perfect  signal  reconstruction  under  both  MC 
and  RIP  criteria.  In  Section  III,  we  first  define  the  probability  of  perfect  recovery.  Next,  we  derive 
the  probability  of  perfect  recovery  in  the  sense  of  both  mentioned  conditions.  Numerical  examples 
are  presented  in  Section  IV. 

Notation:  Throughout  this  paper  bold  uppercase  letters  stand  for  matrices  and  bold  lowercase 
letters  stand  for  column  vectors.  f3[>jds  the  statistical  expectation  and  I  denote  an  identity  matrix. 
In  addition,  ( A)NxM  denote  a  matrix  with  N  rows  and  M  columns  and  {  }>(c  represents  the  real 
part  of  a  complex  quantity.  Further,  ^  ^  denote  the  l2  norm  of  a  and  ||  x||  represents  the  absolute 
value.  Moreover,  <  ~>f  x>  introduce  the  inner  product,  ( >f7  denote  the  transpose  operator,  and  ( '^u 
stands  for  the  Hermitian  transpose. 


II.  CS-Based  MIMO  Radar 


A.  Signal  Model 


In  this  section,  we  describe  the  signal  model  for  our  MIMO  radar  system.  We  assume  that  there 
are  Mt  transmitters,  Mr  receivers,  and  K  targets.  Further,  we  assume  that  each  of  the  targets  con¬ 
tains  Q  individual  isotropic  scatterers.  Let  rjpq  denote  the  transmitted  signal  from  l Ah  transmitter, 
reflected  by  the  gth  scatter  on  the  fc-th  target,  and  received  by  the  p-th  receiver.  Such  a  signal  can 
be  written  as: 


Tlp(Xk))  exp(27T  1  fc(t 


Tip Ofc)))l  exp(2vr 


Tlp(Xk))), 

(1) 


where  wi  ( >f  denote  the  transmitted  waveform,  E  represents  the  average  energy  of  the  transmitted 
signal,  and  fc  is  the  carrier  frequency.  Moreover,  r/p(  >f  and  f['n(  9j.  are  delay  and  Doppler  shift 
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corresponding  to  the  fcth  target,  respectively.  Note  that  {  }wi(t  rip(xk))  exp  2tt  lfct\  stands 
for  the  passband  signal,  (3%  denote  the  reflectivity  of  the  g-th  scatter  on  the  k- th  target,  and  xk  is 
defined  as  the  coordinator  of  the  k- th  target  on  the  range  axis. 

Equation  (1)  can  be  written  as: 


Tip(xk))ex p(27T  1  fc(t  Tip(x fc))) I  exp(27T  l/£(f  Tlp(xk))), 

(2) 


where  Sq  =  j  ®=1 is  the  total  reflectivity  of  the  k- th  target.  Considering  the  narrow-band 
assumption,  we  ignore  the  time  delay  and  rewrite  (2)  as  [12],  [13]: 

K  Mt 

rP{t)  =n  \  Jf' n  {  }®'WexP(2^  1/c(^  rlp(xk)))\  exp(2vr  1  /£(t  TIp(a;fc))) 

fc=l  V  *  z=i 


(3) 


The  baseband  representation  for  fp{t)  is: 

K  r-g-  Mt 

rp(t)  =E[  \  JT3qYI  w*WexP(  2tt  lfcTip(xk))  exp(27r  l/,p(t  rip(xfc))).  (4) 
k= i  V  *  1=1 

Next,  we  transform  (4)  to  discrete  domain  as: 

K  /-gr  Mt 

rP(n)  =n  Jjy^qYI  wi(nTs)ew{  2vr  lfcTiP(xk))  exp(2vr  1  fip{nTs  Tip[xk)))  (5) 

k= 1  V  *  i=l 

We  discretize  the  range-Doppler  space  on  a  M  ±  L  grid  and  rewrite  (5)  in  the  CS  framework  as: 


rp(n)  ~  (cn)lxMi  (A-ri) MtxML^) MLxli 


(6) 


where  c„  =  [wi(nTs),  w2(nTs),  xxx,  wMt(nTs)],  s  = 
elements  of  An  can  be  expressed  as  (al±)n  =  er27ry/  lfrTs 


fI3h  xxx  JWfL  and  the 

'tP(iMi)e2?rVrl^pKs-rMip(a;jtfi))  fQr 


i  =  1,  xxx,  Mt,  k  —  1, 2,  xxx,  ML. 
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Equation  (6)  can  be  expressed  as: 


Applying  the  measurement  matrix,  we  obtain: 

y(n)  =  <f>m( rp)n  =  <pm[ 0,  xxx,  Cn,  >o<x,  0]7VxMt(An)MtxiWL(s)MLxi,  (8) 

where  4>m  (m  —  1, . . . ,  Mm)  denotes  the  m-th  row  of  the  matrix  <X>.  We  can  rewrite  (8)  in  matrix 
form  as: 


y  =  $rp  =  $^(s)Mixi,  (9) 

where  <E  is  the  measurement  matrix  of  size  Mm  ±  N,  —  (C) NxNMt  (A) NMtxML  represents  the 
basis  matrix,  ©  =  'E'E  denotes  the  sensing  matrix  of  size  Mm  ±  ML,  and  s  is  a  A-sparsc  vector. 


B.  Conditions  on  Mutual  Coherence 


MC  of  sensing  matrix  ©  is  defined  as: 


/i(©)  =  max  ^<R6''6"  >  l 


(10) 


where  0,  and  0j  are  different  columns  of  matrix  ©.  We  aim  to  find  the  probability  of  perfect 
recovery  in  terms  of  //(©). 

If  the  following  condition  is  satisfied 

K  <0.5  11 


,(0)L’  (11) 

then  it  is  guaranteed  that  the  OMP  algorithm  obtains  the  true  s.  As  stated  in  [16],  (1 1)  is  a  sufficient 
and  necessary  condition  for  the  sparse  signal  recovery. 
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C.  Conditions  on  RIP 


Suppose  that  for  every  K -sparse  vector  vk,  the  constant  <5  exists  such  that 


(!  &)Jk  j  J  ~  (1  +  5)jk  J- 


(12) 


Let  the  isometry  constant  5  k  is  defined  as  the  minimum  of  all  5  values  that  satisfying  (12).  As 
proved  in  [15],  if  RIP  property  is  hold  for  5k  <  then  it  is  guaranteed  that  the  OMP  algorithm 
converges  to  the  true  s.  Note  that  RIP  is  the  sufficient  condition  for  perfect  signal  recovery. 


III.  Perfect  Recovery  Probability 
A.  Recovery  Probability  in  MC  sense 

Defination:  The  probability  of  perfect  recovery  in  the  MC  sense  is  defined  as: 

a’<05)1  +  mW 

Based  on  the  definition  of  \P,  we  have: 

Mt 

ipuv  =  J\Wi(uTs)ex p(  27T  lfcTlq(xv))  exp(27T  1  fiq(uTs  Tiq{xv ))),  (14) 

1=1 

where  xv  and  ff  are  the  range  and  Doppler  points,  respectively.  As  stated  in  (15),  MC  can  be 

rewriten  as  the  maximum  of  off-diagonal  elements  of  the  Gram  matrix,  i.e.  //  ((-■)  )  =  niax^j  }|[/tJ||  . 

In  the  context  of  CS,  the  basis  matrix  is  typically  pre-determined,  while  the  measurement  matrix 
is  random.  Therefore,  the  behavior  of  MC  depends  on  $  matrix.  In  general,  elements  of  $  are 
independently  chosen  from  the  random  sequence  “}<S>1|  ”  or  a  Gaussian  distribution.  As  we  will 
explain,  applying  Gaussian  distribution  leads  to  high  complexity. 

Let  us  define  the  Gram  matrix  G  as  follows: 

G  =  =  [gi,g2,  xxx,gAr] ,  (15) 
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where  gj  =  [gXj  yx> <gi:j  xxxgNj]T  for  j  =  1,2,  xxx,  N.  Here,  gtj  =  J  k^kj(J  ?(Vv  J  p<t>pq<f>Pk))- 

To  derive  the  probability  of  perfect  recovery,  we  need  to  consider  only  the  upper-right  or  lower- left 
half  of  the  off-diagonal  elements  of  G.  Note  that  g^’ s  are  correlated  RVs. 

Here,  we  assume  that  g0fj  contains  the  upper-right  half  of  the  off-diagonal  elements  of  the  Gram 
matrix.  So,  we  have: 


(16) 


The  goal  of  our  analysis  is  to  quantify  the  probability  of  max^j}  |  1 1  .  For  any  i  §£  j,  |^|fs  are  de¬ 

pendent  RVs.  To  overcome  this  issue,  we  replace  goff  with  the  vector  g  =  \(j\  g2  (j:>,  ...  gN(N-i)/2]T , 
where  g  e  N(/tt,  Cg).  Note  that  2  =  (Jf  ^~1)/2  \\g^f  =  J  ^~1)/2  \fa\f.  Thus,  applying 
eigenvalue  decomposition,  we  have: 


c  = 


' 

A  0 

> 

o 

- 

J 

0  0 

uo 

where 


and 


Thus,  we  have 


U  Un  V 


H 


u 


u  H 

uo 


U: 


H 


u 


\1H 

uo 


u  Un  V=I. 


tci 

1 - 

A  0 

u* 

"I  / 

A  0 

> 

O 

L  u  Uo  v= 

\1H 

u0 

0  0 

uo 

J 

0  0  1 

Finally,  after  whitening  process,  we  have: 


(17) 


(18) 


(19) 


(20) 


g  = 


H 


u 


uo 


(21) 
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Note  that  ^ ^  =  (f  .  |^||)2  =  J  i  It  is  important  to  link  the  probability  density  function 
(PDF)  of  with  that  one  of  maxj}|^|p|  . 

Lemma  1  Let  us  assume  that  the  vector  g  of  size  nR  ±  1  is  given  by  h  =  Kc h  +  h/,  where 
hj  G  Nc(0, R).  Considering  the  fact  that  i£[(h  Ffch)(h  KchH)\  =  E[hfhJ]  =  R, 


we  have  h  G  Nc(  Kch,bR).  Let  Ai,  xxx,  Ar  be  the  nonzero  eigenvalues  of  R  which  means 


M|h||2(s)  — 


1  J 


s  X, 


exp 


sA,iCe||iifh|P 

1  s\i 


Ar+i  =  xxx=  \nR  =  0.  The  moment  generating  function  (MGF)  of  is  obtained  as: 


(22) 


where  the  eigenvector  rq  corresponds  to  the  eigenvalue  A,  (/'  =  1,2,  xxx,  nR). 
Proof:  See  the  Appendix.  ■ 

Now,  equation  (16)  can  be  rewritten  as: 


—  goff  2  ~  max}  lb 


N 


Q  ’  ~  liiCllXV 

V  *7^ 


Hk\ 


/“V 


(23) 


Note  that  the  MGF  of  RV  X  can  be  formulated  as  the  Laplace  transform  of  fx(u),  i.e. 


Thus,  the  PDF  of  X  = 


Mx(s)  —  fx(x)ex p(  sx)dx. 

o 

can  be  found  through  the  following  equation: 


(24) 


fx(x)  = 


C+V- 1° 


exp  (sx)Mx(s)ds. 


27T  1  c— v/— lo_ 

Lemma  2  The  probability  of  perfect  recovery  PMC  is  given  by 


(25) 


PMC  =  Pr(^(0)  >  Vth)  = 

An 


v3.X.L) 


ds1  (26) 


1  c— loo  S  J  "  J  ^ 
where  rjth  —  1/(2  K  1)  denote  the  perfect  recovery  threshold  //(©),  and  the  MGF  Mz{f  is  given 
in  Lemma  1.  ■ 


Proof:  Let  us  define  Y  =  X  /rj  —  ^Jf  lb*  IP/77,  where  77  =  N(N  l)/2.  Then,  the  PDF  of  Y 
is  given  by: 


fr(y )  = 


v2f*(y2 

J_ 

2  y 


v2y 


n 


c+y/— lcxD 


1  c—  \J  —  lo 


exp(sri2y2)Mx(s)ds. 


(27) 
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Replacing  I/77  by  a  curve  fitting  parameter  a,  results  in: 


-yMC 


1 

7TQ!2 

1 

1 

7TQ;2 

T 

1 

■  Vth 

-  C+V— loo 

\ 

sy2  ' 

_  9 

y  exp 

0 

c— V- loo 

J 

or 

'  C+V—  loo 

■  nth 

\ 

si/2  ’ 

2/ exp 

C — y/ —  loo 

0 

/ 

az 

c+C=T°o  Mfs 

)\ 

\  sVth 

-  |  exp 

Mx(s)dsdy 
Mx(s)dyds 
ds. 


(28) 


2vr  1  c-^Too  s  J  "  / 

This  completes  the  proof.  ■ 

We  can  use  (26)  to  determine  the  perfect  recovery  probability.  Note  that  only  one  integration  is 
involved  in  the  obtained  solution. 


B.  Recovery  Probability  in  RIP  sense 

In  this  subsection,  we  first  define  the  probability  of  perfect  recovery  and  then  derive  it  as  a 
function  of  5  and 


/V 


Defination  The  probability  of  perfect  recovery  in  RIP  sense  is  defined  as: 

PRIP(S,  Jkj)  =  Pr(l  S  yK  \  ~  U  ~  1  +  5  yK  |), 


vv 


VV 


(29) 


where  5  is  a  constant,  vfc  denote  any  A'-sparse  vector,  and  U  =  K ■ 

Here,  we  define  #  =  [rif)1,  xxx, and  $  =  (f>p,  xxx,  <fiTMrJT ,  where  4>i  (i 

1,  2,  >oo<,  Mm)  denote  the  i-th  row  of  matrix#. 

Based  on  (12),  we  rewrite  the  middle  term  of  inequality  (14)  as: 


©V^  \  =  ##VA' 


(30) 


/  7 

where  is  a  A'-sparse  vector,  and  n  =  # # v K  =  (ri\.  xxx,  nr .  'XXXKMrn)r.  In  addition,  n,  = 
f  pJi  Vp^p.  As  a  result,  we  have  ^ )vK  =  f  |^>jC|p,  where  c  stands  for  the  constant 


vector  c  =  f  ^rJ[  Considering  the  independency  property  and  recalling  the  center  limited 

theorem,  we  conclude  that  is  a  Gaussian  RV.  In  what  follows,  we  find  the  mean  and 

variance  of 
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Since  different  elements  of  matrix  $  are  independent  and  identically  distributed  Gaussian  RVs 


with  equal  variance  a2,  the  mean  value  of  can  be  expressed  as: 


E[V 

The  variance  is  obtained  as: 


=  E[]J  \  fac\f}  =  c H  jn  77, 


[ 


c. 


(31) 


Var[^K2]  =  E 


Mr, 


ri 


=  E 

M, 


{E[j>vK*])2 

]'^(</>,c)2  +  2r/(0<c)l  (E[^>wK2})2 

m 

n  Bd^cif]  v« jo2 + evm>  <  c)]. 


(32) 


In  deriving  (32),  the  following  relation  is  used:  )  f  ;=rJl  ||0MmcIP  [  =  r/(0,  c)2  +  2 r/(0  <  c),  where 


v(<i>,  c) = n  ii^cif 


(33) 


2—1 


and 


C2 

'-'A 


V(<t><c)  =  II  I^CIN^C|P- 


(34) 


k= 1 


The  PDF  of  is  fully  characterized  by  equations  (31)  and  (32). 


Lemma  3  The  probability  of  perfect  recovery  PRlp(fi, 


p  (Vv)=  g 


la, 


7 


2 +  2) 


/V 


r + q  i 


is  given  by: 


/i(2  Vv’ 


r. 


(35) 


where  a2  denote  the  variance  of  RV  U,  5  is  a  given  constant  in  RIP  definition,  ()( >)  is  the  Gaussian 
Q-function,  and  ^  stands  for  an  arbitrary  A'- sparse  vector.  In  deriving  (36),  we  used  the 
relation  erfi(x)  =  erfc(  lx)/  1,  where  erfc  =  2  ^jj3xp(  t2)/  ndt.  ■ 
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Proof:  The  probability  PRIP(S:  obtains  when  matrix  ©  satisfies  the  defined  RIP  condi¬ 
tion  in  (12)  for  given  vector  vk  and  5.  Considering  the  PDF  of  we  have: 

r  (i+<5)IIvkH! 

'(1-<5)I|vx||1 

1  1 


°Vv)  = 


fu(u)du 


_  exp  )  ^  ^ 

2ttou  2  V  “  ‘/  4fff 


=  Q 


k  2  +  2) 


/a. 


V 


r + q  ^ 


erfi 

1 

at, 


1  ,  P 

u  + 


2  *1  2al 

(2  $Jk  J)  \  ■ 


2m2 


P-+5)\\vk\\1 


■i— 5)l|vj<riii 


(36) 


This  completes  the  proof. 


IV.  Numerical  Results 

In  this  section,  we  present  numerical  results.  The  targets  are  assumed  to  fall  on  the  grid  points. 
Throughout  our  simulations,  the  carrier  frequency  is  assumed  to  be  1  GHz,  Mt  =  Mm  =  10,  and 
K  /  }1,  3,  5|  .  In  addition,  the  grid  size  is  set  to  105  ±  91  and  10  compressed  measurements  is 
used  at  each  receiver. 

Figure  1  shows  the  probability  of  perfect  recovery  in  MC  sense.  Here,  the  elements  of  matrix 
$  are  drawn  from  a  zero-mean  Gaussian  distribution.  According  to  equation  (11),  radar  is  able 
to  successfully  detect  K  =  6  targets  with  the  probability  of  90%,  while  it  can  successfully  detect 
K  =  5  targets  with  the  probability  of  99.9%.  In  addition,  decreasing  the  variance  value  leads  to 
better  performance.  For  example,  changing  cr2  from  0.7  to  0.1,  improves  the  recovery  probability 
from  90%  to  96%. 

In  Fig.  2,  we  present  the  probability  of  perfect  recovery  in  RIP  sense.  Here,  we  consider 
three  different  values  for  K .  It  is  observed  that  increasing  K  improves  the  probability  of  perfect 
recovery.  Also,  it  is  seen  that  increasing  the  probability  of  perfect  recovery  leads  to  increase  in  5 
value. 
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Appendix 


Proof  of  Lemma  1 

We  first  apply  a  whitening  process  for  the  RV  g.  Let  g  denote  the  outcome  of  the  whitening 
process.  Then,  the  mean  and  variance  values  of  g  can  be  written  as: 

uf  Kg 

■m  = 1 


n  =  E[  g]  = 


U 


H 


TJ  H 

u0 


(37) 


uf  Kg 


and 


Var(g)  =  E[(g  /x)(g  i*)H]  = 


A  0 

0  0 


I  ’ 


(38) 


where  A  is  a  diagonal  matrix.  Further,  /i  =  (/x1?  /ij,  ynR),  where  /i,  =  K uf  g.  As  a  result,  we 
have  g  E  Nc  (//,  Var(g)).  Note  that  all  of  the  elements  of  g  are  independent  RVs  with  distribution 
hi  E  Nc(fii,  Aj).  Furthermore,  we  have  ^ ^  =  f  ri=1  where  ||r),|fs  are  Ricean  RVs 

with  the  following  PDF 


t  t  \  2  Xi  ^  K+IKIP 

f\~9i\{Xi )  =  -J7exP  J  - ^ - 


2  IK  If* 

A* 


(39) 


Here,  /0(>)  denote  the  modified  Bessel  function.  So,  the  MGF  of  h  2  is  found  as: 

r, 


B[exp  ,s-nifelfl]  =  S[/  exp 
/  i=  1  J=  1 


■IfelfD 


XXX 


00  ,  ,  /•  2 r 2xi  \  K  +  IKIP 

ax  i  >oo uxr  /  exp  sx{  [  — —  exp  I  - - - 

o  /  ,  A*  /  A* 


exp 


#=i 

(1  sK)xi  +  IK  If 


2|KIK 

A,- 


Aj 


2  IK  IK 


Aj 


dxt 


(40) 


2y< 


'Ll 


1  sAj  q  A 


exp 


1  sAj 


exp 

■sAj|KlP 

1  sAj 


y2  +  IK  IP 

Aj 


2  Vi  |  K 


Aj  1  sAj 


dyi 


This  completes  the  proof. 


12 


References 


[1]  E.  Candes  and  T.  Tao,  “Decoding  by  linear  programming,”  IEEE  Trans.  Inf  Theory ,  vol.  51,  pp.  4203-4215,  Dec.  2005. 

[2]  J.  A.  Tropp  and  A.  C.  Gilbert.  “Signal  recovery  from  random  measurements  via  orthogonal  matching  pursuit,”  IEEE  Trans. 
Inf.  Theory,  vol.  53,  pp.  4655-4666,  Dec.  2007. 

[3]  S.  S.  Chen,  D.  L.  Donoho,  and  M.  A.  Saunders,  “Atomic  decomposition  by  basis  pursuit,”  Soc.  Indus.  Appl.  Math.  Rev.,  vol. 
43,  pp.  129-159,  Mar.  2001. 

[4]  R.  G.  Baraniuk,  “Compressive  sensing,”  IEEE  Signal  Process.  Mag.,  vol.  24,  pp.  118-124,  July  2007. 

[5]  E.  J.  Candes,  M.  B.  Wakin,  “An  introduction  to  compressive  sampling,”  IEEE  Signal  Process.  Mag.,  vol.  25,  pp.  21-30,  Feb. 
2008. 

[6]  J.  Wang,  S.  Kwon,  and  B.  Shim,  Near  optimal  bound  of  orthogonal  matching  pursuit  using  restricted  isometric  constant, 
EURASIP  J.  Adv.  Signal  Process.,  vol.  2012,  no.  1,  pp.  1-8,  2012. 

[7]  M.  A.  Herman  and  T.  Strohmer,  “High-resolution  radar  via  compressed  sensing,”  IEEE  Trans.  Signal  Process.,  vol.  57.  pp. 
2275-2284,  June  2009. 

[8]  Y.  Zhang,  “Theory  of  compressive  sensing  via  Zi-minimization:  A  non-RIP  analysis  and  extensions,”  CAAM  Technical  Report, 
Rice  University,  TR08-11,  Houston,  TX,  2008. 

[9]  J.  Wang  and  B.  Shim,  “On  the  recovery  limit  of  sparse  signals  using  orthogonal  matching  pursuit,”  IEEE  Trans.  Signal 
Process.,  vol.  60,  pp.  4973-4976,  Sept.  2012. 

[10]  M.  A.  Davenport  and  M.  B.  Wakin,  “Analysis  of  orthogonal  matching  pursuit  using  the  restricted  isometry  property,”  IEEE 
Trans.  Inf.  Theory,  vol.  56,  pp.  4395-4401,  Sept.  2010. 

[11]  S.  Satpathi,  R.  L.  Das,  M.  Chakraborty,  “Improving  the  bound  on  the  RIP  constant  in  generalized  orthogonal  matching 
pursuit,”  IEEE  Signal  Process.  Lett.,  vol.20,  pp.  1074-1077,  Nov.  2013. 

[12]  Y.  Yao,  A.  P.  Petropulu,  H.  V.  Poor,  “Measurement  matrix  design  for  compressive  sensing-based  MIMO  radar,”  IEEE  Trans. 
Signal  Process.,  vol. 59,  pp.  5338-5352,  Nov.  2011. 

[13]  S.  Gogineni  and  A.  Nehorai,  “Target  estimation  using  sparse  modeling  for  distributed  MIMO  radar,”  IEEE  Trans.  Signal 
Process.,  vol.59,  pp.  5315-5325,  Nov.  2011. 

[14]  A.  P.  Chandrakasan,  V.  M.  Stojanovic,  “Design  and  analysis  of  a  hardware-efficient  compressed  sensing  architecture  for  data 
compression  in  wireless  sensors,”  IEEE  J.  Solid-State  Circuits,  vol.  47,  pp.  744-756,  Mar.  2012. 

[15]  M.  A.  Davenport  and  M.  B.  Wakin,  “Analysis  of  orthogonal  matching  pursuit  using  the  restricted  isometry  property,”  IEEE 
Trans.  Inf.  Theory,  vol.  56,  pp.  4395-4401,  Sept.  2010. 

[16]  J.  Wang  and  B.  Shim.  ”A  simple  proof  of  the  mutual  incoherence  condition  for  orthogonal  matching  pursuit”  submitted  for 
publication. 

[17]  I.  S.  Gradshteyn  and  I.  M.  Ryzhik,  Table  of  Integrals,  Series,  and  Products,  6th  ed.  San  Diego:  Academic  Press,  2000. 

[18]  V.  K.  Rohatgi,  An  Introduction  to  Probability  Theory  Mathematical  Statistics.  Wiley:  New  York,  1976. 


13 


Pr(n(0)<n 


^th 


Fig.  1 .  Probability  of  perfect  recovery  in  MC  sense  with  different  values  of  variance  cr^ 
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Fig.  2.  Probability  of  perfect  recovery  in  RIP  sense  for  given  J-f-sparse  scenarios,  where  K  £  {1,  3,  5}. 
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